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Based on the analogy with non-minimal SU(2) symmetric Wu-Yang monopole with regular metric, the solution 
describing a non-minimal U{1) symmetric Dirac monopole is obtained. In order to take into account the curvature 
I coupling of gravitational and electromagnetic fields, we reconstruct the effective metrics of two types: the so-called 

associated and optical metrics. The optical metrics display explicitly that the effect of birefringence induced by 
^ ' curvature takes place in the vicinity of the non-minimal Dirac monopole; these optical metrics are studied analytically 

, and numerically. 

0^ 
(N . 

1. Introduction 
O 

The Dirac monopole as a specific static spherically symmetric solution to the minimal Einstein-Maxwell equations 
^ has became a subject of discussion in tens papers, reviews and books (see, e.g., [TJ [51 |2l IH 13 El Ej ) ■ The motion of 
{^jTji massive and massless particles, which possess electric charge or are uncharged, is studied in detail (see, e.g., [319] 
'— '■ and references therein). In the paper [TD] we introduced and discussed the SU{2) symmetric Wu-Yang monopole of 
, ' the new type, namely, the non-minimal monopole with regular metric. Since the non-minimal Wu-Yang monopole 
^ is effectively Abelian, it is naturally to consider the corresponding analog of that solution in the framework of non- 
f — minimal electrodynamics. Mention that non-minimal models with magnetic charge have been discussed earlier (see, 
• e.g., iJJiliSj), but the exact analytical regular solution obtained here as direct reduction to the C/(l) symmetry is 
^ ■ the new one. 

. The second novelty of the presented paper is the investigation of photon dynamics in the vicinity of the Dirac 

' monopole accounting the curvature coupling of the gravitational and electromagnetic fields. In the presence of 

1 I non-minimal interaction (induced by curvature) the master equations for electromagnetic and gravitational fields 

in vacuum can be rewritten as the master equations in some effective anisotropic (quasi)medium [131 114j . This 
means that two effective (optical) metrics can be introduced |15[ I16[ I17j , so that the photon propagation in vacuum 
interacting with curvature is equivalent to the photon motion in the effective space-time with the first or second 
optical metric, depending on the photon polarization. Even if the real space-time has a regular metric, the optical 
^ metrics can be singular, admitting the interpretation in terms of the so-called "trapped surfaces" and "inaccessible 
zones" [TH1[TH]. We discuss this problem in Sect. 3. Numerical modeling of the photon orbits, presented in Sect. 4, 
supplement our conclusions. 

2. Master equations and background fields 

We consider the LagrangiarH 

+ ^F^kF'" + i7^^'=™"F,fcF„„ (1) 

OTT Z Z 

to describe the background gravitational and magnetic fields in the framework of non-minimal Einstein-Maxwell 
model with the non-minimal susceptibility tensor 



^[9 9 -9 9 ) - J^^-K +4(/i: .g + H g -Kg - H g ) 



(2) 



linear in curvature. As usual, i?*'^"™ is the Riemann tensor, R""^ is the Ricci tensor, R is the Ricci scalar, Fmn is 
the Maxwell tensor. This non- minimal susceptibility tensor can be obtained from the general one (see [Sj), when the 
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coupling constants gi , (72 , 93 are chosen as follows qi ~ —q < 0, q2 — 4:q, qs — —6q. The ansatz for the space-time 
metric is 

ds^o) = N{r)dt^ - ~ {dO^ + sin^ 0dip^) . (3) 

Mention that the choice gu9rr = ^ 1 in our ansatz is supported by the result, obtained in [10] for the given relations 
between the coupling constants qi , (72 , qs- The equations of electrodynamics [15[ I20j 

WkH''' = , VfcF*'*^ = (4) 

with constitutive equations 

Tjik /^ikmn -r^ /^ikrnn ^ / ^irri kn ^in km\ , ^ikmn /r-\ 

are associated with the Lagrangian where F**'"' is the dual Maxwell tensor and H^'' is the induction tensor. The 
equations ^ with (O are satisfied identically, when the potential of static spherically symmetric electromagnetic 
field out of the point-like magnetic charge fj, is of the form 

Ak^S'^A^^-St ^1(1- cos 9). (6) 

The corresponding strength field tensor has only one non-vanishing component 

Fe^ = -^sinO , (7) 

which does not depend on the non-minimal coupling parameter q. Thus, the well-known solution with the magnetic 
field of the monopole type 



= F*'''Uk ^ Si.^ , B{r) = ^^B^, = ^ , (8) 



satisfies the non-minimal Maxwell equations ([U, ([5]), The equations for the gravity field 

R^k - Ir 9^k = STTT^f , (9) 



obtained by direct variation procedure from the Lagrangian ([T|) are non-minimally extended, since the effective 
stress-energy tensor T^''^^' has the form 



T^f^ -T^~l\T^k'~ ^k' + (^nk'l ■ (10) 
The quantities , xj;^'^ , T^^^ and T^j^'^ are given by 



^ik — ^9ik^mn-f --tim-fk i U-Lj 

Tlk^ = R - \R^kFmnF^'^ - V VKi^™„i^"") + ^ V, Vfc (F,„„i^"") , (12) 
Tlk^ = -\9^k [V™V,(^^™"^^'„) - i?,™F""i^'„] - F''\RaFkn + RklF^n) - R"''^FmFkn 

- ^V'V,(F,„i^fe") + [WdFknF'") + WkiFnF'")] , (13) 

tI^^ = ^g^kR^^'^'^FmnFis - ^F'«(F, "i?fc„i, + F.'^Rms) - ^V™V„(F, "F,™ + F^^F, ") . (14) 

We are happy to stress that with the ansatz ([7]) for the metric ([3]) the system of equations ^ transforms into one 
equation 

-^'(l + S+^(l-3v)-l-^-3V- (15) 
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whose solution is 



N{r) = 1 + 



r2(K - 4Afr) 
2(r4 + Kq) 



(16) 



Here k is the convenient constant, k — 87r/i^ . This solution is regular in the center (-/V(O) = 1) and satisfies the 
asymptotic condition A^(oo) = 1 . If the mass M is less than its critical value M((.j.it) , where 



A 



Q \ J 2 

the metric ^ with (fT6|) has no horizons as in [TO] . 



. 48 g ^ 
1 + — ^ + 1 



(17) 



3. Electrodynamic description of the photon propagation 

Linear electrodynamics allows us to consider the dynamics of test photons in terms of microscopic field strength fa- 
and induction hik , which satisfy the equations 

Wkh''' = 0, = C*'=™"/™„ , Vfc/"^ = 0. (18) 

Clearly, the microscopic tensor fmn is the analog of the macroscopic Maxwell tensor F,„„ , describing the background 
electromagnetic field, and h^'' is the analog of the tensor iJ'*^ (compare ([4]), ^ and (fT8|) ). We have to stress that 
the tensor of material coefficients C'*'''™" in the formula (fT8|) describing the photon propagation, is the same one 
as in ([5|) for the background field. Thus, the non-minimal interaction influences the photon by two ways: first, via 
the constitutive equations (see ([5]), second, via the gravitational field, whose potentials (fTB|) depends on the 

parameter of non-minimal coupling q . 

3.1. Associated metrics 

The tensor C'*'^™" can be represented algebraically as a quadratic combination of the so-called associated metrics 
q^^f") [1511121 



(19) 



where G(^a){0) and ft compose a set of parameters, introduced in [T7]. When the medium is spatially isotropic, this 
decomposition requires one associated metric; in the case of uniaxial spatial symmetry one needs two associated 
metrics [161 117) . In the framework of our model vacuum interacting with curvature can be regarded as an uniaxial 
quasi-medium, where the radial direction is the selected one. Indeed, let us use the standard definitions of the 
dielectric permittivity tensor e*'^ , magnetic impermeability tensor {n~^)pq , and tensor of magnetoelectric coeflacients 
i^p", given by [UdD] 

e — Z(Vfe(7„0 , )pq——-ripik^ Vmnq , — r]pikC Un- 



Here U = S^/vN is the velocity four-vector, associated with the magnetic charge (at rest), r]pik = epikmU'^ and 
Cpifcm is the Levi-Civita tensor. In the spherically symmetric model under consideration the tensor is equal to 
zero identically, other tensors have the following non-vanishing components: 

e|l^£; = 2C'■Vt = l + 27^'^Vt, (21) 
e^^el^el^l + 2n\, = 1 + 27^^;, , (22) 
(A.-l); = l + 27^^^. , (23) 
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For the metric ^ with we obtain 



et 



7^' 



(r4 + a4)3 



o 4 7 

3a r 



(25) 



7?'"* 



(r4 + a4) 



4 10 



(26) 



7^' 



(r4 + a4)3 



[12gMr5 - 5aV* - 20gAfaV + 3a^] , 



(27) 



A new parameter a with dimensionaHty of length is defined as foUows: — nq. Since £g = and (// ^)9 = (/^ ^)S^, 
the medium can be considered as uniaxial. Moreover, due to the equality "R-^^^g = T^^^gt the relation £_lijl± = 1 takes 
place. The tensor (7'*""" can be reconstructed algebraically as 



^ 2/1 



gim(A)^kn{A) _ ^m(A)^fcm(A) 



■ 7 



gim{A)_girn(B) \ lkn{A) _ gkn(B) 



^[A)_ 



')(» 
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where the associated metrics g^^'^^'f and are 
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Mil 



Ai± Aim , 



(28) 



(29) 



Here X^^^ = d'^y/N{r) is the radial tetrad space-like four- vector and Ajij 
representation (|28p is a particular case of the decomposition (fT9|) with 



(30) 

U^Um is the projector. The 



- = '^ll 



i - 1 _ ^ll^-L 
7 e-LMll 



(31) 



Mention that (jr''^(^) can be obtained from g^^'^^'^ by a formal replacement of the symbol e by the symbol /i. 
3.2. Optical metrics 

In the approximation of geometrical optics ^2lj the potential four-vector Ai and the field strength tensor can be 
represented as 



respectively, where $ is the phase, a/ is a slowly varying amplitude, and is a wave four-vector: 
km = . 

In the leading-order approximation the Maxwell equations reduce to 
C kikfYiO-n — . 



(32) 
(33) 
(34) 



This system of linear equations with respect to a„ admits non-trivial solutions, when four components of ki satisfy the 
Fresnel equation (see, e.g., [15] for details), which is usually called the dispersion equation. Using the decomposition 
((28l) with ((29)) and l|30)) in (p4|) after some algebra we obtain the dispersion equations for two different principal cases. 

(i). The first case relates to the propagation of electromagnetic wave with non- vanishing longitudinal polarization, 
i.e., with a,. 7^ 0. The Fresnel equation admits the solution 



^ + krk'' + ^—{kgk^ + k^k'^) ^ , 



(35) 
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where to = U^ki is a frequency. This dispersion relation can be reconstructed as an eikonal-type equation 

5""(^)fc.A;™ = 0, (36) 

where the first associated metric (P^ . 

(ii). The second case is characterized by = 0. The Fresnel equation admits the solution 

^2 1 

-^ + krk'' + (fcefc^ + fc^fc'^) = 0, (37) 

£±/i|| 

which can be rewritten as 

5™(s)fc,fc„ = 0, (38) 



where is the second associated metric ((30)) . 

When the electromagnetic wave propagates along the radial direction, i.e., when kg — k^p — and fc^fc^ = —k"^, 
the dispersion relations ((35| and ((37|) give the same result, = fc^c^ , providing the phase velocity to be equal to 
speed of light in standard vacuum. The velocity of propagation in the transversal directions (i.e., when kr = 0) 
depends on the polarization: when 7^ the phase velocity of the wave is V(i) = ^/f^, when = it is 

V, 



(2) 

4. Photon dynamics 

The associated metrics obtained in Subsection 3.1 are the optical ones, i.e., they satisfy the eikonal equation. This 
means that the motion of a photon influenced by the non-minimal interaction can be described by null geodesic lines 
in the effective space-times with optical metrics (jr""^"*) or , depending on polarization. Three metrics: the 

metric of space-time ([3]), the first and second optical metrics (|29|) and pO|) . can be written in general form 

dsfs) = N{r)dt^ - - r^Y^E){r){de^ + sin^ Od^^) , (39) 

where the index E takes the values 0, A, B and F(o)('') = l- For the sake of simplicity we consider below the 
monopole with M = 0. It is characterized by three interesting features: first, N{r) > 1 for arbitrary r, thus, N{r) 
does not reach zero; second, N(r) reaches its maximum at r = a, the maximal value being N(^„iax) 
third, A^(0) — 1 and curvature invariants are regular in the center r = . For this particular model we have 

i+9^-7e+e ' ^ ^ 

where ^ is dimensionless positive quantity ^ = Kqjr^ . The polynomial 1 — 11^-1-37^^-1-^'^ has no real positive zeros, 
i.e., y(^)(r) ^ 0. The equalities Y(^)(r) = 00 and simultaneously y(_B)(r) = take place, when 1-1-9^ — 7^^-1-^3 = 0, 
i.e., at C = 6 ~ 1-85 and C = 6 ~ 5.25. Finally, Yi^B)(r) = 00, when C = 6 ~ 0.19 and C = ^4 ~ 0.54. Thus, 
despite the metric (is^Q^ is regular everywhere, the first and second optical metrics have singular points, when the 
angular functions ^(a) ('') ^^nd ^(b) i^r) vanish. Physical interpretation of such (dynamic) singularities will be done 
in the next Subsection in the course of description of radial, tangent and inclined particle motion. 

4.1. Photon trajectories 

Let us consider null geodesic lines in the effective space-times with optical metrics and (|30p . The geodesic 
equations have the standard form 

(P'x^ dx^ dx^ „ 

^+^^V)^^=0. (42) 
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Figure 1: Rays for optical A-metric at the impact parameter J > J^^j^ ^ 0.794659: (a) J = 2, (b) J = 1, (c) J = 0.8, (d) 
J = 0.79466. 



where the ChristofFel symbols ^jn^E) : = 0' ^ ^-^^ constructed on the base of the corresponding metric. As in 
classical case (see, e.g., [55^) the photon moves in the plane, and one can choose this plane as the equatorial one, 
9 = Tr/2. Other integrals have the form 

dip J dt 1 

d^ 



r^YrE){r) ' dr N{r) ' 



dr\ 
d^l 



= 1- 



J^N{r) 



(43) 
(44) 

dr dr ^ ^ ^ 

where the quantity J is the impact parameter. As a consequence of and the following relation takes place: 



providing the relation 



9(E)r. 



dr^ = , 



1 



—r'Y^E)ir)-N{r) 



(46) 



The formulas are the non-minimally modified formulas for the well-known integrals of motion 22j (they 

coincide with them, when Y(^E){f) — 1 only). Let us emphasize seven interesting details of these modified integrals 
of motion. 

(I) Radial motion is non-sensitive to the non- minimal coupling. Indeed, when 6 = ^ and ip = const, the 
parameter J vanishes and the relation ^ = ±1 is valid for each E = 0, A, B . In this case the velocity of photon 
coincides with speed of light in standard vacuum. 
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(II) When the motion is not pure radial (inchned), the condition = describes the extrema of the curve r{ip) , 
giving the distance of the closest approach. When the curvature interaction is absent, i.e., q — 0, we obtain from 
()46p the standard condition for the minimal distance tq: — j'^N{ro) — 0. In the non-minimal model Y(yi) and 
Y(B) a-re the functions of radius (see (HOI) and (dH)), and the set of solutions to the equation Yi^e){'>''^Y{e) ~ J'^N) — 
is much more sophisticated than in the minimal model. 

(III) Formally speaking, the optical metrics become non-Lorentzian, i.e., are of the signature H 1 — h, when 

Yi^E){f) is negative, and are singular at the points, where Y(£;)(r) = 0. From the dynamic point of view the negative 
values of Y(£;)(r) are, nevertheless, admissible. Indeed, the angular integral of motion in (I43p remains consistent, but 
the direction of rotation becomes opposite; the right-hand-sides of the equations (f44|) and (I46|) remain positive, as it 
should be. 

(IV) The functions Y(£;)(r) change the sign at the points, where Y(^E){f) = or Y(^E){f) = °o ■ The first possibility 
can be realized for the B-ray only, when ~ nq/r'^ takes the values 1.85 or 5.25. In this case the equation (|46p 
gives the condition ^ = 0, indicating the point of extremum for the curve r{{p) (see item III). The conditions 
Y[E){f) = can be realized at ^ = 1.85 and ^ = 5.25 for the A-ray, and at ^ = 0.19 and ^ — 0.54 for the B-ray. 
At these points the angular frequency ^ vanishes (see (|43p ) and the radial velocity ^ coincides with speed of light 
in the standard vacuum (see (|44l) ). Since the relation ^ — co holds at these points, the rays are directed radially. 

(V) The relation ^ = is self-contradictory, when Y(^e) is negative (see (011)). This means that tangent waves 
can not propagate, and we can indicate this zone as inaccessible one for them. For the tangent A-ray, there is only 
one inaccessible zone at 1.85 < ^ < 5.25, for the tangent B-ray the second one appears at 0.19 < ^ < 0.54. The 
corresponding boundary spheres can be indicated as trapped surfaces (see, e.g., [HI [19] for details). 

(VI) For A- and B-rays there are specific "critical" values of the impact parameters 

-^crit'^^' = max ) . 
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Figure 3: Rays for optical B-metric at the impact parameter J > j'^f^ ~ 5.29385: (a) J = 6, (b) J = 5.4, (c) J = 5.294, 
and (d) J KJi^l, J = 5.29. 



If J > Jcrit, the corresponding rays do not reach the surface with Y(B)(r) = oo. Mention that in minimal case such 
"critical" values of the impact parameters do not exist. 

(VII) Since the space-time metric of the non-minimal monopole is regular, i.e., N{r) ^ and N{r) 7^ oo, in the 
t — r cross-section there is no singularities. 



4.2. Numerical modeling of the photon trajectories 

Pictures presented in the Fig. 1-4 illustrate the dependence r{'^) for rays, which are null geodesies of the optical 
metrics g^^^ and 5^^^ . The rays start from infinity, the impact parameter being equal to the height of the starting 

point of curve. We put k = 1, g = 10. In the Fig. 1, 2 two circles relate lo r — Ri = ^ Kq/1.85 and r = 
i?2 = ^K,q/5.25; in the Fig. 3, 4 three circles relate to r = R3 = ^Kq/0.19, r — R4 = ^ Kg/0.54, and r = Ri, 
respectively. The picture (a) in the Fig. 2 contains one self-intersection point, the picture (d) in the Fig. 3 and the 
picture (a) in the Fig. 4 contain two self-intersection points. The region inside the circle r = R\ is inaccessible zone 
for the inclined B-rays. These figures show that the non- minimal Dirac monopole can be considered as a scattering 
center for inclined rays. Scattering laws for A- and B-rays are different confirming the birefringent character of 
non-minimal interaction. 
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Figure 4: Rays for optical B-metric at tlie impact parameter J < J^^^^ « 5.29385: (a) J = 5, (b) J = 3, (c) J = 1, (d) 
J = 0.1. 



5. Conclusions 

Qualitative and numerical analysis of the photon orbits in the vicinity of non-minimal Dirac monopole with regular 
metric has demonstrated the following interesting features. 

1. Propagation of the electromagnetic waves in the vicinity of non-minimal Dirac monopole is characterized by 
birefringence, induced by curvature, i.e., the phase velocities of waves depend on their polarization. Two different 
optical metrics should be introduced to describe two principal states of polarization. 

2. The metric of the non-minimal Dirac monopole, obtained and discussed in this paper, is regular, thus, all the 
singularities of the optical metrics have a dynamic origin and are supported by the non-minimal (curvature induced) 
interaction of the gravitational and electromagnetic fields. 

3. The points of self-intersection, the points of the closest approach, the reverse points, etc. in the photon 
trajectories can be recognized and catalogued for different combinations of the values of the impact parameter and 
the parameter of non-minimal coupling. We discussed here only the principal pictures. 
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